DESY 98-195 
KA-TP-20-1998 
tiep-ph/9812298| 
October 2001 (rev.) 

Top Dipole Form Factors and Loop— induced 
CP violation in Super symmetry 

W. Hollik a ' 6 , J.I. Illana c , S. Rigolin d 
C. Schappacher a , D. St6ckinger a 

a Institut fur Theoretische Physik, Universitat Karlsruhe, 
D-76128 Karlsruhe, Germany 

b Theoretical Physics Division, CERN, 
CH-1211 Geneva 23, Switzerland 

c Deutsches Elektronen-Synchrotron DESY, 
D-15738 Zeuthen, Germany 

d Departamento de Fisica Teorica, Universidad Autonoma de Madrid, 
Cantoblanco, E-28049 Madrid, Spain 

Abstract 

The one— loop Minimal Supersymmetric Standard Model (MSSM) contributions to 
the weak and electromagnetic dipole form factors of the top quark are presented. Far 
from the Z peak, they are not sufficient to account for all the new physics effects. In 
the context of the calculation of the process e + e~ — > tt to one loop in the MSSM, we 
compare the impact on the phenomenology of the CP-violating dipole form factors of 
the top quark with the contribution from CP-violating box graphs. Some exemplifica- 
tive observables are analyzed and the relevance of both the contributions is pointed 
out. The one-loop expressions for the electromagnetic and weak dipole form factors 
in a general renormalizable theory and the SM and MSSM couplings and conventions 
are also given. 



1 Introduction 



The investigation of the electric and magnetic dipole moments of fermions provides deep 
insight in particle theory. The measurement of the intrinsic magnetic dipole moment (MDM) 
of the electron proved the correctness of the hypothesis of half-integer-spin particles |J and 
is one of the most spectacular achievements of quantum field theory predictions. More precise 
studies of electron and muon showed afterwards the presence of anomalous contributions to 
the MDM (AMDM). Their predictions constitute one of the most spectacular achievements 
of Quantum Field Theory and imply very accurate tests of the quantum structure of the 
Standard Model (SM). The measurements of the (g e — 2) and (g^ — 2) available |2| are in 
perfect agreement with the SM predictions to several orders in the perturbative expansion of 
the theory (cf. and references therein). Furthermore, with the expected precision at the 
E821 Brookhaven experiment |3J it will be possible to improve the previous measurement of 
(g M — 2) by a factor 20. Therefore the MDMs can be used, together with the precision tests 
at the Z resonance from LEP and SLC and the new results of LEP2 and TEVATRON, to 
set bounds on possible new physics effects beyond the SM |J. 

The importance of the analysis of the electric dipole moment (EDM) of elementary and 
composite particles is intimately related to the CP violating character of the theory. In 
the electroweak SM there is only one possible source of CP violation, the Sqkm phase of 
the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix for quarks ||. Currently the only 
place where CP violation has been measured, the neutral K system, fixes the value of 
this phase but does not constitute itself a test for the origin of CP violation [^]. On the 
other hand, if the baryon asymmetry of the universe has been dynamically generated, CP 
must be violated. The SM cannot account for the size of the observed asymmetry ||. In 
extended models (beyond the SM) many other possible explanations of CP violation can 
be given. In particular, in supersymmetric (SUSY) models || [H| CP violation can appear 
assuming complex soft-SUSY-breaking terms. Two physical phases remain in the GUT 
constrained MSSM [|H], [H], ET2J, enough to provide the correct size of baryon asymmetry in 
some range of parameters |TI|. But the most significant effect of the CP violating phases 
in the phenomenology is their contribution to the EDMs |H| . Unlike the SM, where the 
contribution to the EDM of fermions arises beyond two loops [pL5 1 , the MSSM can give a 
contribution already at the one-loop level |J. 



The measurements of the neutron, electron and muon EDMs fll6| , [17[ constrain the phases 
and the supersymmetric spectrum in a way that may demand fine tuning (supersymmetric 



CP problem): either the SUSY particles very heavy (several TeV [[L8[]) or the phases are 
of 0(1CU 2 ) ||. Very large soft-SUSY-breaking masses are unappealing as it seems natural 
to demand the SUSY spectrum to be at the electroweak scale.[] On the other side, the ex- 
perimental constraints can be met taking general universal soft-SUSY-breaking terms and 

1 Moreover if the SUSY spectrum is in the TeV region this could also give rise to relic densities unacceptably 
large. 
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vanishing SUSY CP phases. In this case the CP violation is originated via the usual SM 
CKM mechanism and the super symmetric spectrum affects the observables only through 
radiative contributions. In this scenario it is difficult to contruct models in which the SUSY 



phases naturally vanish [191 and at the same time provide some other non-standard mech- 
anism for explaining electroweak baryogenesis. There exist also ways of naturally obtaining 
small non-zero SUSY CP phases which leave sufficient CP violation for baryogenesis [ |12|| . 
But, in general, they do not lead to observable differences from the SM. In models with 
potentially observable predictions one has to relax the assumption of soft-term universality. 
Several attempts have been made, following this direction, to use CP violation from top- 
squark mixing: a complex parameter A t would yield large CP violating effects in collider 

)6|.0 Besides, due to renormalization-group relations, 

One can satisfy the experi- 



processes involving top quarks |20 



the phase of A t is constrained by the EDM of the neutron [22 



mental constraints due to cancellation among the different components of the neutron EDM 
(constituent quarks and gluons), the SUSY phases can still be kept of 0(1) and the SUSY 
spectrum at the electroweak scale satisfying the experimental bounds ||23|| . In |24|, it is 
shown that large CP violating phases in the MSSM are compatible with the bounds on the 
electron and neutron EDMs as well as with the cosmological relic densities. In view of all 
these arguments we keep our analysis completely general and consider the SUSY CP-phases 
as free parameters. 

Some attention has also been payed to the study of possible CP violating effects in the 
context of R-parity violating models f26| . In this class of models new interactions appear 
providing extra sources of CP violation (still preventing fast proton decay). They can explain 
the CP violation in the K system (with no need of the CKM phase) without introducing 
anomalous Flavor Changing Neutral Current (FCNC) contributions p7 |. In the following 
we restrict our discussion to the simpler case of models with conserved R-parity. 

In addition a certain amount of investigation has been devoted to the analysis of weak 
dipole moments (WDM). The WDMs are defined, in analogy to the usual DMs, taking the 
corresponding on-shell chirality-flipping form factors of the Zff effective vertex. The SM 
one-loop contribution to the anomalous weak magnetic dipole moment (AWMDM) has been 
calculated for the r lepton and the b quark in [^, |2"9f . The CP-violating weak electric dipole 
moment (WEDM) is in the SM a tiny three-loop effect. The WDMs are gauge invariant 
and can directly connected to physical observables. While for the r case, using appropriate 
observables |28| , |30| , 3jJ , an experimental analysis is feasible, for the b case the situation is 
complicated by hadronization effects ||32|| . The SM predictions are far below the sensitivity 



reachable at LEP J28| but non-standard interactions can enhance these expectations (2HDM 
33[, MSSM g) especially for the (CP-violating) WEDMs (2HDM |5j, leptoquark models 



3§, MSSM 0]). The experimental detection of non-zero AWMDM or WEDM of heavy 



fermions, at the current sensitivity, would be a clear evidence of new physics beyond the SM. 



Large non-SM CP violating top-quark couplings could be probed at high energy colliders like the NLC 
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In the perspective of the next generation of linear e + e~ colliders, we extend the previous 
analyses on the WDMs to consider the t quark dipole form factors (DFF) and, in particular, 
the CP-violating ones. In |38| an independent analysis of the t quark EDFF and WEDFF 
can be found. Since the t quark is very heavy one expects this fermion to be the best 
candidate to have larger DFFs. Beyond the Z peak (s > Amf) other effects are expected 
to give contributions to the physical observables. In fact the DFFs in a general model are 
not guaranteed to be gauge independent. An exception is the one-loop MSSM contribution 
to the CP-violating electromagnetic and weak DFFs (EDFF and WEDFF) which do not 
involve any gauge boson in internal lines. Anyway, although this contribution is indeed 
gauge invariant any CP-odd observable will be sensitive to not only the DFFs but to the 
complete set of one-loop CP-violating diagrams involved. In this work we present the full 
calculation of the expectation value of a set of these observables in the context of the MSSM 
and compare the size of the different contributions. 

The paper is organized as follows. In Section 2 we present the general effective vertex 
describing the interaction of on-shell fermions with a neutral vector boson. The definitions 
and the generic expressions of the DFFs for all the contributing topologies at the one- 
loop level are given. In Section 3 we briefly discuss the cancellation of the dipoles in the 
supersymmetric limit. The numerical results for the t quark EDFFs and WEDFFs at a/s = 
500 GeV are presented in Section 4. In Section 5 we evaluate specific CP-odd observables 
for the t quark pair production in e + e~ colliders to one loop in the MSSM and compare the 
influence of the t EDFF and WEDFF with that of the CP-violating box diagrams. Our 
conclusions are presented in Section 6. In the Appendices one can find the definition of the 
one-loop 3-point tensor integrals employed in the expressions for the dipoles as well as the 
SM and MSSM couplings and conventions that have been used. 



2 The dipole form factors 
2.1 The Vff effective vertex 

The most general effective Lagrangian describing the interaction of a neutral vector boson 
V with two fermions can be written, using at most dimension five operators, as a function 
of ten independent terms: 



C 



Vff 



V{x)^{x) 



7 M (flv - QaIb) + i {9m + i#E75) 



+i d v Ots + h7TP75) 



ty(x) 



+ (\d»V v {x))m(x) 



9fiu (igs + gpls) + o-^ (kjtm + #TE75) 



$?(x). 
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Table 1: C, P, T properties of the operators in the effective Lagrangian of Eq. (|7p. Their 
chirality flipping behavior is also displayed. 



Operator 


Coefficient 


P 


CP 


T 


Chirality Flip 






+ 


+ 


+ 


NO 






— 


+ 


+ 


NO 


<— > 




+ 


+ 


+ 


YES 












YES 




fl'TS 


+ 






YES 




#TP 




+ 


+ 


YES 


(9 ■ V) 


#S 


+ 






YES 


(id ■ V) ^75^ 


#P 




+ 


+ 


YES 


( W) tt^tf 


#TM 


+ 


+ 


+ 


YES 




fi'TE 








YES 



The first two coefficients, i.e. gy and #a, are the usual vector and axial-vector couplings. 
They are connected to chirality conserving dimension four operators. All the other coef- 
ficients in Eq. (|l|) multiply chirality flipping dimension five operators and can receive a 
contribution only through radiative corrections in a renormalizable theory. The operators 
associated to gy, g&, gu, 9p, 9tm and #tp are even under a CP transformation. The presence 
of non vanishing g#, g$, (?te and g^s yields a contribution to CP-violating observables. In 
Table |1] we summarize the C, P, T and chirality properties of each operator introduced in 
the effective Lagrangian. 

By Fourier transform of Eq. (|1|) one obtains the most general Lorentz structure for the 
vertex Vff in the momentum space: 



rVff 



7 M (/v - /a7s) + {q ~ q)p Urn + i/fiTs) + P» 0/s + /p7b) 

+(q - qfa^ (/ TS + i/TP7s) + p"cv (1/tm + /te7s) , (2) 

where q and q are the outgoing momenta of the fermions and p = (q + q) is the total incoming 
momentum of the neutral boson V. The form factors f\ are functions of the kinematical 
invariants. Actually they are more general than the coefficients g±. In fact any operator of 
dimension higher than five added to the Lagrangian of Eq. ([!]) is related to a new coefficient 
But every new coefficient can contribute only to the ten independent form factors /j 
introduced in Eq. (|2]). The parameters g^ and the form factors can be complex in general. 
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Their real parts account for dispersive effects (CPT-even) whereas their imaginary parts are 
related to absorptive contributions. 

It is possible to lower the number of independent form factors in Eq. (0) by imposing on- 
shell conditions on the fermionic and/or bosonic fields. For instance, in the case of on-shell 
fermions, making use of the Gordon identities: 

2m f u^v = jw (q - q~Y v + u i (q + q) I/ a fJ,u v\, 

2m f uYl5 v = {u (q + qYl5 v + u i (q - 9)^0^75 u}, 

= [u{q + qY v + ui{q-q) v a^v], 

= {u{q-qY 1 c,v + ui(q + q) u a^^ 5 v} , (3) 

one can eliminate /tm, /te, /ts an d /tp from the effective Lagrangian. The number of 
relevant form factors can be further reduced taking also the boson V on its mass shell. In 
this case the condition p M e M = automatically cancels all the contributions coming from 
fs and fp. The same situation occurs for off-shell vector boson V when, in the process 
e + e _ — ► V* — > //, the electron mass is neglected. Therefore fs and fp will be ignored in 
the following. With all these assumptions the Vff effective vertex for on-shell fermions is 
conventionally written as: 



T v "{s) = ie 7m \VY{S) - 4» 75 1 + <v(9 + ?) 



1 o " 75 

2m j e 



(4) 



where e and m/ are respectively the electric unit charge and the mass of the external fermion. 
The form factors in Eq. (Q) depend only on s. As mentioned above Vj (s) and Ajf (s) param- 
eterize the vector and axial- vector current. They are connected to the chirality conserving 
CP-even sector. The form factors aX{s) and dj(s) are known respectively as anomalous 
magnetic dipole form factor (AMDFF) and electric dipole form factor (EDFF). They are 
both connected to chirality flipping operators. In a renormalizable theory they can receive 
contributions exclusively by quantum corrections. The EDFFs contribute to the CP-odd 
sector and constitute a source of CP-violation. 

The dipole moments (DM) are defined taking the corresponding vector bosons on shell, 
s = My. For V — 7 one gets the usual definitions of the photon anomalous magnetic dipole 
moment (AMDM) and electric dipole moment (EDM). The definitions of electric charge, 



magnetic and electric dipole moments consistent with our convention for the covariant 
derivative ( |B2"D are respectively: 



charge = -e V?{0) = e Q f , (5) 
— 1 



MDM = (y7(0) + oj(0)), (6) 

2m f J 



EDM = cT f {0) . (7) 

Thus the AMDM of a fermion is aj(0) = —Qf{g/ — 2)/2 with gf being the gyromagnetic 
ratio. The axial- vector coupling Aj vanishes. For V = Z, the quantities af = (M§) 
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Class I: Class II: Class III: 




Figure 1: The one-loop Vff diagrams with general couplings. 

and dj = g?J(M§) define the anomalous weak-magnetic dipole moment (AWMDM) and the 
weak-electric dipole moment (WEDM). 3 



2.2 One— loop generic expressions of the dipole form factors 

All the possible one-loop contributions to the a^(s) and d^(s) form factors can be classified 
in terms of the six classes of triangle diagrams depicted in Fig. [I|. The vertices involved are 
labeled by generic couplings for vector bosons — W^, WL fermions and 

scalar bosons according to the following interaction Lagrangian: 

£ = ieJ(Wl u W^V - W^WlVy + V^W^Wu) + ieG jk V»<S>) H M $ fc 
+{e* f {S jk - PjklB^k^ + eKjkVVWQj + h.c.} 

if 75)*; • (8) 



Every class of diagrams is calculated analytically and expressed in terms of the couplings 
introduced in (§D and the one-loop 3-point integrals C$ (see App. [A]). The result is given in 
the 't Hooft-Feynman gauge. 



[Class I]: vector-boson exchange: 



2m f 47r L 7k~i 



(I) - TZ { 4m f E Re MP ( V ff V fk* + AfjA®*" 



3 Massless and neutral fermions may have magnetic moments. The usual parameterization in (|j) must 
be generalized by the replacement e aj(s)/2mf — > /i (s). The dipole moments of neutrinos are then given 
by n v (My). 
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2C^~ — 3C-l~ + Cq 



-^W^-^)][2Cf-C ] WI } 



(9) 



(I) 



47T 



(Ot/(0*> 



kji 



jki 

+A%\V?}V?> + Af s A%)] [2C}- - Cf 
-4i: mJ m[^>(V^4'>--4] l /<f) 

-^(^V/f - ^g*)] [2Ct - C ] J (10) 

In the case of gluon exchange one has to substitute a for a s and for T\ (the SU(3) 
generators) and take A^- = 0. The sum over the index / yields a color factor Cp = 4/3. 
The gluon contribution to the CP-violating form factor c?Jf(I) vanishes in general. 

[Class II]: fermion exchange and two internal vector bosons: 
-(H) 



4(f). 
2m 4 



a 
An 



{2m f £ Re [j(V®vjt» + A$Af>)] [AC+ + C+ 
jki 



kji 



-6Z™MAvjPvjt»-A$Af* 



jki 



1 ikjii 



(II) = 



a 
An 



{2m f ZH J ( V fl A T + A flrf )m ) 
jki 

+6j2milm\j{v}pA ( f f 



ACl 



kji 



A$V<t>) 



jki 



1 Ikji) 



(11) 



(12) 



• [Class III]: scalar exchange: 
a){s) 



2m f 



(HI) = 



djf( f 



"(HI) 



-2 J2 m k Re [vjP (5 y S? fc - P^) 



-^(^ - ^)] [^ + + cr]^} 



An 



{-2m f Y: im [y£ v) (p^; fe + s l3 p; k ) 

jki 

+A ( j k \Sij S* k + Pij Pi* k ) 



2(7^ — C l 



kji 



+2j2m k lm[v} k V) (P lj S: k - S^) 

jki 



(13) 
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[Class IV]: fermion exchange and two internal scalars: 



ct + cr] } 
1 1 ikiii 



4(f), 

2m i 



OL ( r 

(IV) = --{2m / ^Re[G ifc (^*+P i ^) 



2Ci 



jkl 



+ ^m,Re[G ifc (^ - PflPSa)] [^t ~ C ]J (15) 



jki 



kji 
Co 



(14) 



kji 



4(, 



-(IV) 



a 



47T L 



{2m / X;i m fG ifc (^+P i ,5^ 



2C+ 



jM 



- J2 mM[G jk (SjiPki - PjiS* kl ) 
jki 



kji 



(16) 



[Class V+VI]: fermion exchange, one vector and one scalar internal boson: 



2m/ 
4( S ) 



(V + VI) 
(V + VI) 



a 

Air' 



2 £ Re K jk {Vjt ] Sl + A$P$ C+ + C\ 



(k) 



jki 



kji 



a 
4n 



2 J2 ImlK^fp; + A%%)\ [Ct + CI 



(*), 



(17) 
(18) 



In Eqs. (PHlSl) the shorthand notation [C]kji stands for the 3-point tensor integrals 
C(—q,q,Mk,Mj,Mi). The integrals appearing in the previous Eqs. are UV and IR fi- 
nite. All the expressions are proportional to some positive power of a fermion mass, either 
internal or external, consistently with the chirality flipping character of the dipole moments.^ 
For class V and VI diagrams the mass dependence is hidden in the product of the Yukawa 
couplings Sij (Pij) and the dimensionful parameter Kij. Hence, the heaviest fermions are 
the most promising candidates to have larger DFFs. Eqs. (|9|- [T8|) also show that, in general, 
the DFFs for massless fermions are not vanishing but proportional to masses of fermions 
running in the loop. The SM cancellation of the massless neutrino DFFs is only due to the 
absence of right-handed neutrinos. Finally, notice that all the contributions to the EDFFs 
are proportional to the imaginary part of certain combinations of couplings. A theory with 
real couplings has manifestly vanishing EDFFs. 



3 Cancellation of the dipole moments in the supersym- 
metric limit 

A general Vff interaction is restricted to the form @ by Lorentz invariance. Since the 
Lorentz algebra is a subalgebra of supersymmetry, this interaction is even more constrained 

4 This is as expected when applying the mass-insertion method: to induce a flip in the fermion chirality 
one introduces a mass in either one of the internal fermion lines, picking a mass term from the propagator, 
or in the external fermion lines, using the equations of motion for the free fermion. 
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in a theory with unbroken supersymmetry. In supersymmetric sum rules are derived that 
relate the electric and magnetic multipole moments of any irreducible N = 1 supermultiplet. 
Applied to the Vff interaction between a vector boson coupling to a conserved current and 
the fermionic component of a chiral multiplet, these sum rules force the gyromagnetic ratio 
to be gf = 2 and forbid an electric dipole moment: 

a) = d) = 0. (19) 

The Lagrangian of the MSSM is supersymmetric when the soft-breaking terms are removed. 
For non zero value of the fi parameter the Higgs potential has only a trivial minimum. 
Therefore to keep the particles massive and supersymmetry preserved at the same time the 
choice fi = is necessary. Then the Higgs potential is positive semi-definite and it has 
degenerate minima corresponding to v = v i = v% (tan/3 = 1). The value of Ma = follows 
from such a configuration. Finally the value of the common v is fixed by the muon decay 
constant. 

In this supersymmetric limit the above mentioned sum rules are valid and the magnetic 
and electric dipole form factors have to cancel. To verify our expressions we checked this for 
the AWMDM of the b quark [El. Choosing the parameters A b = A t = M 2 = M 3 = 0, /i = 0, 



tan P = 1 and Ma = 0, the SM gauge boson contribution to a™ |29| is indeed cancelled by 
the MSSM correction including the two Higgs doublets: the gluon and gluino contribution 
cancel among themselves and the neutralinos and charginos cancel the gauge boson and 
Higgs contributions. A similar check has been performed for the electric and weak-electric 
dipole form factors of the t quark. 



4 SM and MSSM predictions for the top quark dipole 
form factors 



4.1 SM 

The electroweak contributions to the magnetic and weak-magnetic dipole form factors for 
off-shell gauge bosons are gauge dependent. The pinch technique |41| could be used to con- 
struct gauge-parameter independent magnetic dipoles in the class of gauges [[4^, |33| but 
the prescription is not unique and these quantities cannot be observable by themselves. 
The QCD contributions (gluon exchange) to the t (W)MDFF are gauge independent and 
comparable in size to the electroweak predictions at yfs = 500 GeV in the 't Hooft-Feynman 
gauge due to the large mass of the t quark. 

There is no contribution to the electric and weak-electric dipole form factors to one loop 
in the SM. 
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4.2 MSSM 



The triangle diagrams with SUSY particles are gauge independent by themselves (no gauge 
or Goldstone bosons involved) but the ones including Higgs scalars and gauge bosons in the 
loop are not sufficient to keep the gauge invariance in the case of the magnetic dipole form 
factor. The CP-violating dipole form factors (electric and weak-electric) for which the Higgs 
sector of the MSSM is irrelevant, on the other hand, can be regarded as gauge independent 
quantities. These form factors have also been considered in [38] and the numerical results in 
are in agreement with the ones presented here. Since the MSSM contains a CP-conserving 
Higgs sectoif] and the SM provides contributions to the (W)EDMs beyond two loops, the 
relevant diagrams are the genuine SUSY graphs of classes III and IV. 

We investigate the SUSY contributions to the t electric and weak-electric dipole form 
factors performing a parameter scan for which a fixed value a/s = 500 GeV has been chosen.0 
It is important to point out that the region of the supersymmetric parameter space that 
provides the maximal contributions varies with s due to threshold effects. We scan the mass 
parameters M 2 and in a broad range and the CP-violating phases: ip^, V 9 * an d <Pi (see 
Figs. ^] and |3|). The gluino mass (jrig = M 3 ) is given by the GUT constraint (fTq). We 
adopt a fixed value for the common scalar quark mass rriq = 200 GeV: this is a plausible 
intermediate value; larger values decrease the effects. Finally, the moduli of the off-diagonal 
terms in the t and b mass matrices are also chosen at fixed values = = 200 

GeV, to reduce the number of free parameters. The results are expressed in t magnetons 
jj, t = e/2m t = 5.64 x 10~ 17 ecm. 

The MSSM contributions to the top quark EDFF and WEDFF are analyzed below. 



4.2.1 Neutralinos and t scalar quarks 

They provide typically small contributions but quite sensitive to the value of both the phases 
involved, ip^ and (pi (Fig. §). 

The results are larger for low tan/3 since the chirality flipping mass terms are domi- 
nated by the t quark, yielding a term proportional to m t cot (3. A term proportional to the 
neutralino masses is also present as well as a negligible one proportional to mfetan/3. The 
contributing diagrams belong to the classes III and IV for the Z case and only to class IV 
for the 7 case, as the neutralinos do not couple to photons. 

As a reference we take the representative values M 2 = = 200 GeV and tp^ = = 
7r/2. For these inputs the results are 

d][x°] = ( 0.080 + 0.081 i) x 10~ 3 [i t , (20) 
df[x°] = (-0.324 + 0.223 i) x 10~ 3 m . (21) 



5 A one-loop non vanishing contribution to the WEDM is possible in a general 2HDM |35 . 

6 We use the running coupling constants evaluated at y/s — 500 GeV, a s = 0.092, a = 1/126. 
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ReK[x°]} 



Im{^[x ]} 




Re{df[x ]} 




/ 




Im{df[x ]} 



2tt 



I / +0.2 



2tt 




Figure 2: Neutralino contribution to the real and imaginary parts of the t EDFF and WEFF 
[in lCT 3 /! t units] in the plane ipi — y? M for tan/3 = 1.6 and the reference values M 2 = = 
| m *. fl | = 200 GeV aty/s — 500 GeV. 
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4.2.2 Charginos and b scalar quarks 

As in the case of the neutralinos, the influence of chargino diagrams is enhanced for low 
tan j3. 

The results depend very little on (p~ b and mainly on ip^, the maxima being close to 
<Pn = ±7r/2. For increasing tan/3 the dependence on pi~ h grows, as it comes with a factor 
proportional to m b tan (3. 

The chargino contributions are the most important ones. In Fig. ^| the dependence on 
Im(/i) and M 2 is displayed for tan/3 = 1.6. The only relevant CP-violating phase here has 
been set to the most favorable case, ip^ = tt/2 (the negative values of Im(/i) correspond to 
tpp = —tt/2). The symmetry with respect to Im(/i) = in Fig. |3| reflects the (almost) inde- 
pendence of ifi, here set to vr/2.|] The same does not happen for the neutralino contributions 
for a fixed value of pi ^ 0, tc in the plane M 2 — lm(/j,). The plots of Fig. ||| exhibit a tendency 
to decoupling of the supersymmetric effects for increasing values of the mass parameters 
f44| . The isolines for a couple of masses of the lightest charginos and neutralinos in the 



same plane are also given for orientation. The current LEP2 experimental lower limits are 
m^o > 25 GeV and m-± > 90 GeV @. 

The chargino contributions for the values M 2 = — 200 GeV and tp^ — 7r/2 are 



= (0.869 - 1.870 i) x 10~ 3 fi t , (22) 
df [^±] = (0.793 - 2.524 i) x 10~ 3 [i t . (23) 

4.2.3 Gluinos and t scalar quarks 

Their effect is roughly proportional to |tt4^| sin p- t times a chirality flipping fermion mass, 
either m t or M 3 (the gluino mass). It is damped for heavy gluinos circulating in the loop 
and also for large scalar quark masses due to decoupling. Both terms have opposite sign to 
lm(m t LR ) and the one proportional to the gluino mass dominates. 

The result for M 2 = 200 GeV and p- t = -tt/2 is 



d][g] = (0.457 + 0.170 i) x 10~ 3 /i* , (24) 
df[g] = (0.155 + 0.059 i) x 10" 3 [i t . (25) 

7 All the contributions flip sign when the set ((p^, fit Vb) 1S rotated by it. They vanish accordingly when 
all the phases are zero. See Fig. || for illustration. 



13 



Re{d7[x±]} 



M, 



500 



400 



300 




ImK^]} 



M, 



500 



400 



300 



200 -F^ 



M, 



500 



400 



300 



200 



100 



500 


-250 


250 500 






Im(^) 




Re{df[x±]} 







M, 



500 



400 



300 



200 



100 



500 


-250 


250 500 






Im(^) 




Im{df[x±]} 





-500 -250 250 500 

Im(^) 




-500 -250 



250 500 
Im(^) 



Figure 3: Chargino contribution to the real and imaginary parts of the t EDFF and WEFF 
[in lCT 3 ^ units] in the plane M 2 — Im(yu) for tan/3 = 1.6, \m b LR \ = 200 GeV and ip~ b = n/2 
at y/s = 500 GeV. The lower (upper) solid isolines correspond to m~± = 90 (200) GeV and 
the dashed isolines to m^o = 25 (100) GeV. 
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4.2.4 Total contribution 

In view of these results, we establish a set of SUSY parameters and phases for which nearly 
all the contributions sum up constructively at yfs = 500 GeV. Our choice is 

Reference Set #1 : tan/3 = 1.6 

M 2 = \n\ = rriq = |m* LR | = \m\ R \ = 200 GeV 

W = ~ t Pt = ~ ^l^i (26) 

for which the t EDFF and WEFF reach the values: 

dj = (1.407 - 1.618 i) x 10" 3 fi t (27) 
df = (0.624 - 2.242 i) x 10~ 3 fx t , (28) 

at t/s = 500 GeV. To illustrate how the maximum effects appear we display in Fig. || the 
individual contributions as a function of \fs. The masses of the supersymmetric partners in 
the loops are such that there are threshold enhancements in the vicinity of y/s = 500 GeV. 



5 CP— odd observables 

One needs to go beyond the Z peak to produce t quark pairs in e + e~ collisions and this implies 
to account for several features. The photon-exchange diagram is no longer suppressed in 
this regime and therefore one needs to separate the contributions of electromagnetic and 
weak form factors f45| . Moreover, not only vertex- but also box-diagrams correct the tree 



level process to one loop. The former corrections are parameterized by the electromagnetic 
and the weak vertex form factors but the latter demand the introduction of new form factors 
according to the more general topology of the process. For a realistic theory, one expects 
that any CP-odd observable will depend not only on the CP-violating effects due to vertex 
corrections included in the electric and weak-electric dipole form factors but also on possible 
CP-violating box contributions. This is indeed the case of the supersymmetric models. 

We concentrate on supersymmetric CP violating effects in t-pair production at e + e~ 



colliders. CP conserving MSSM one-loop contributions to e + e~ — > ff are discussed in [46 
To our knowledge, only the electric and the weak-electric form factors have been considered 
so far to parameterize these effects.^ Our purpose is to evaluate the expectation value 
of several observables in the context of the MSSM with complex parameters to one loop. 
Eventually we compare the contribution from the electric and weak-electric form factors to 
these observables with the CP-violation effects coming from the box contributions. 

Our starting point is an initial CP-even eigenstate in the c.m.s. (laboratory frame). This 
is the case for unpolarized e + and e~ beams but only a good approximation in the case of 



A similar analysis for hadron colliders has been recently presented in Rcf. 
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Figure 4: The different contributions to the t (W)EDFF [in 10 3 fi t units] for the reference 
set of SUSY parameters of Eq. Ml 
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longitudinally polarized beams (unless they are ideally 100% polarized): neglecting the 
electron mass, chirality conservation allows only interactions of equal helicity states that 
project into each other under a CP transformation. 

Consider the process e + (p + ) + e~(p_) — > £(k + ,Si) + t(k_,s 2 ) (pair-production of 
polarized t quarks). The decay channels labeled by a and c act as spin analyzers in 
t + t — > a(q + ) + c(q_) + X. The momenta and polarization vectors in the overall c.m.s. 
transform under CP and CPT as follows:^ 



CP : p± -> -p T = p± CPT : p ± -> p T = -p ± 

(29) 



k± — > — k T = k-t k-t — > k T = — k± 



q± -» -q=F q± 

Si <-> s 2 Si <-> — s 2 



From the unit momentum of one of the £ quarks in the c.m.s. (say k + ) and their polar- 
izations (si, s 2 ) a basis of linearly independent CP-odd spin observables can be constructed 
45|. They are classified according to their CPT parity. The spin observables are related 



to more realistic (directly measurable) momentum observables based on the momenta of the 
top decay products P5| . The polarizations can be analyzed through the angular correlation 
of the weak decay products, both in the nonleptonic and in the semileptonic channels: 



t(k+) -> &(q b )^had(qx) , 

t(k+) -> % b )£ + (q+)z/, (£ = e,n,T) 



(30) 
(31) 



and the charged conjugated ones.[^] As m t > Mw + m*, the t quark decays proceed pre- 
dominantly through Wb. Within the SM the angular distribution of the charged lepton is a 
much better spin analyzer of the t quark than that of the b quark or the W boson arising 
from semileptonic or nonleptonic t decays The dimensionless observables are easier to 
measure, for instance the scalar CP-odd observables |H5|: 



A 1 
A 2 



q+ x q 

|q+ x q- 
(q+ + q- 



[CPT-even] 
[CPT-odd] 



(32) 
(33) 



or the CP-odd traceless tensors 451: 



T- 
j — 



Q 



lq+ 

(q+ 



q 



iq+ x q- 



q-Ji^ 



x q_| 
" q-)i 



+ (i <-> j) 



[CPT-even] 
[CPT-odd] 



(34) 
(35) 



9 T means reflection of spins and momenta. 

10 The leading QCD corrections to e + e~ — > tt, that include a gluon emission, have a very small effect on 
the t spin orientation [b0| . 
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The reconstruction of the t frame is not necessary for the momentum observables. The ob- 
servables A 2 and Qij do not involve angular correlations as they could be measured consid- 
ering separate samples of events in the reactions e + e~ — > aX and e + e~ — > aX. Nevertheless 



it is convenient to treat them in an event-by-event basis |45 . 

One may obtain additional CP-odd observables from combinations of the standard ten- 
sors by multiplying them with CP-even scalar weight functions to maximize the sensitivity 
to CP-violating effects (optimal observables). Neglecting quartic terms in the CP-violating 
form factors, A,, the differential cross section can be written as da = d<7 + J2i^A cr i- It 
has been shown that the observables given by Oi = do~{/do~o have maximal sensitivity to 
the CP-violating terms A, |5^| . The CP-odd and CPT-even (odd) observables are sensi- 
tive to the dispersive (absorptive) parts of the CP-violating form factors. These include 
the real (imaginary) parts of the electric and weak-electric dipole form factors as well as 
CP-violating contributions from box topologies. 

Since the t quark is a heavy fermion, the CP conjugate modes tiXh and t^tn are produced 
with a sizeable rate. This allows to construct the following CP-odd asymmetry |53], plj 



This asymmetry is related to the one that can be measured through the energy spectra of 
prompt leptons in the decays t — » W + b — > i^Vfb and its conjugate. The W + is predom- 
inantly longitudinally polarized and, assuming the standard V— A interaction, the b quark 
is preferentially left-handed. The W + is mostly collinear with the t polarization and so is 
the £ + anti-lepton. Above the ti threshold a £ + coming from a tn has more energy, due to 
the Lorentz boost, than one produced in a tz decay. The same happens for the conjugate 
channel, the l~ from a is in average more energetic than the one from the i^. Therefore, 
in the decay of the t^tn the anti-lepton from tn has a higher energy E + , while in the decay of 
the pair tjJz the lepton from II has a higher energy E-. Thus the asymmetry A is sensitive 
to the energy asymmetry of the leptons [Qj or b quarks [^| in the final state. 

We will ignore possible CP violation in the t or i decays. CP-violating t decays in 
supersymmetry have been considered in |56[|. We evaluate the influence on the CP-odd 



observables of the vertex and box diagrams through (ideal) spin observables. The expectation 
value of the realistic momentum observables given above will also be presented (assuming 
SM top decays) for comparison with experimental capabilities. 



6 MSSM full contribution to CP— odd observables 



The one-loop differential cross section for polarized t-pair production in the MSSM [[46 
involves the box diagrams indicated in Fig. |5| besides the vertex graphs of Fig. |l| The class 
with vector boson exchange contains only SM contributions ([eZZt], [i/ e WWb]) whereas the 
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Figure 5: Relevant generic box graphs for the MSSM at one loop. 



Table 2: CP-odd spin observables and the coefficients for the expectation value of their 
integrated version at ^/s = 500 GeV, where only the CP-violating dipole form factors are 
taken into account. 



i 


CPT 


Oi 


a 


b 


Cl 


c 2 


c 3 


c 4 


1 


even 


& 


s 2)y 


N| 


N| 


0.526 





1.517 





2 


even 






N| 


L| 


-0.465 





-0.061 





3 


even 




X s 2;)z 


N| 


T| 


0.708 





0.144 





4 


odd 


to 


~ S 2Jx 


T| 


T| 





0.930 





0.123 


5 


odd 


to 


— $2)2 


L| 


LI 





-1.417 





-0.287 


6 


odd 


(«* 


X S 2 ) y 


L| 


TJ. 





0.263 





0.758 



other one is purely supersymmetric ([ex°X°^]; [&eX a^D an( i provides CP violating effects. 
The box diagrams with Higgs boson exchange are proportional to the electron mass and are 
neglected in the whole calculation (they are CP even in any case). When we refer to one- 
loop calculation in the following, the QED as well as the standard QCD corrections to the 
tree level process are excluded: they need real photonic and gluonic corrections to render an 
infrared finite result and constitute an unnecessary complication as they are CP-even and 
do not affect qualitatively our conclusions. 



6.1 Spin Observables 



A list of CP-odd spin observables classified according to their CPT properties is shown in 
Table |2|. Their expectation values as a function of s and the scattering angle of the t quark 
in the overall cm. frame are given by 



(O)ab 

da 



2da 



E + E 

s^,sj=±a,±b s£,S2=±b,±a 

E da(st,s 2 ). 

±s*,±s* 



da(s* 1 ,s 2 ) O 



(37) 
(38) 
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The directions of polarization of t and t (a, b) are taken normal to the scattering plane (N), 
transversal (T) or longitudinal (L). They can be either parallel (f) or antiparallel (|) to the 
axes defined by z = k + , y = k + x p + /|k + x p + | and x = y x z. Notice that 

(0 2 ) ah = -L £ M<,£) O 2 (39) 

s*,sj=±a,±b 

and the same for the average of any CP-even quantity. If the information of the t scattering 
angle is not available one can consider integrated observables 



E + E 

s*,s*=±a,±b s*,s*=±b,±a 



a(s\,s* 2 )0. (40) 



The contributions to the CP-odd observables are linear in the t EDFF and WEDFF and 
in the CP-violating parts of the one-loop box graphs. The shape of the different dipole 
contributions to these observables is depicted in Fig. |6|. The coefficients of the dipoles in the 
linear expansion of the integrated spin observables, 

(0) ah = ^ ( Cl Re[d]} + c 2 lm[d]} + c 3 Re[df } + c 4 lm[df}) , (41) 

are shown in Table ^| for their integrated version at */s = 500 GeV. They are model inde- 
pendent. 

In Fig. [7] we compare the contributions of dipoles and boxes to the spin observables, for 
two reference set of parameters. The first set was given in fl26D and is the one that has been 
shown to enhance the dipole effects. The second one is: 

Reference Set #2 : tan/3 = 1.6 

M 2 = =rriq = = \m b LR \ = 200 GeV 

<Pn = ft = <fl = V 2 , ( 42 ) 



where, compared to (|2q) , only the CP-violating phases have been changed. The shape of 
the solid and dashed curves is the same in all expected. Their different size is 

due to the contributions to the normalization factors coming from self-energies, A(W)MFFs 
and other CP-even corrections. The plots show that for the Set #1 the MSSM box graphs 
contribute (in the case of CPT-even observables) to CP violation in the process e + e~ — > tt 
by roughly the same amount and with a different profile than the EDFFs and WEDFF of 
the MSSM. For nearly all cases this eventually results in a large reduction of the value of 
any CP-odd observable with respect to the expectations from the dipoles alone, as it will be 
shown below. Other sets of SUSY parameters, that do not enhance the dipole contributions, 
can provide instead smaller dipole form factors but larger observable effects. This is the case 
of the Set #2, that yields much smaller values for the real parts of the dipole form factors 
at = 500 GeV than we had before in Eq. (|28|), 

dl = (0.547- 1.889 i) x 10" 3 fx t , (43) 
df = (0.362 - 2.319 i) x 10~ 3 fx t . (44) 
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Figure 6: Dipole contributions to the expectation values of the spin observables at y/s = 500 
GeV for different values o/(Re[d7], I m M7], Re[df], Im[df]) in m units. 
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Figure 7: Expectation value of the spin observables [in 10 -3 units] for the two reference 
sets of SUSY parameters, assuming for the cross section: (i) tree level plus contribution 
from (W)EDFFs only (solid line); (ii) one loop including all the vertex corrections and the 
self-energies (dashed line); (Hi) complete one loop (dot-dashed line). 
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Table 3: Ratio r = (O) / y (O 2 ) [in 10~ 3 units] of the integrated spin observables at y/s = 
500 GeV for the two reference sets of SUSY parameters. The left column excludes the box 
corrections and the right one comes from the complete one-loop cross section for e + e~ —>■ tt. 



% 


CPT 




a 


b 


Set #1 


Set #2 


1 


even 


(«* 


-s* 2 


)y 


N| 


N| 


1.216 


-0.231 


0.207 


-1.394 


2 


even 




x < 


X 


N| 


L| 


-0.755 


-0.489 


-0.053 


0.318 


3 


even 


( s ; 


x s* 2 


)z 


NT 


T| 


1.184 


0.625 


0.090 


-0.598 


4 


odd 


(s* 


-< 


X 


Tt 


T| 


-1.230 


-1.421 


-1.888 


-2.217 


5 


odd 






)z 


L| 


H 


2.550 


2.739 


3.823 


4.216 


6 


odd 


fa 


X s* 


)y 


L| 


T| 


-1.683 


-1.751 


-2.050 


-2.216 



As a consequence, the CPT-even observables receive their contributions mainly from the 
box diagrams (Fig. |7|) for the Set #2. 

The previous arguments are reflected in Table where the ratio r = {0)/yj (O 2 ) is shown 
for the integrated spin observables. We compare the result when only the self energies and 
vertex corrections are included (left column) with the complete one-loop calculation (right 
column). The shape of the observables as a function of the t polar angle (Fig. 0) is such 
that signal for CP-violation r is still not very large for a couple of CPT-even observables 
with the Set #2 but it is much better for all the others. This illustrates that the dipole form 
factors of the t quark are not sufficient to parameterize observable CP-violating effects and 
the predictions can be wrong by far. 

As a final remark, note that (O5) = A, the asymmetry defined in Eq. ([36]). Attending to 
Table it happens to be the most sensitive observable to the imaginary part of the EDFF 
and also the best one to test CP violation for our choice of SUSY parameters (Table |3J). 
The observables (sj — s 2 ) can still give information on CP violation when the polarization 
of one of the t quarks is not analyzed.^ The sensitivity of the single-spin polarization to 
CP-violation is indeed worse: ((s* - s* 2 ) x ^ z ) ah = 2((s*)^) a for a(b)= N|( j),TT(I),LT(I). 

6.2 Momentum Observables 

Consider now the decay channels labeled by a and c acting as spin analyzers in t + i — ► 
a(q + ) + c(q_) + X. The expectation value of a CP-odd observable is given by the average 

11 Anyway a comparison between two samples, one with polarized t and the other with polarized f, is 
necessary to build the genuine CP-odd observable of Eq. ([38]). 
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Table 4: Ratio r [in 10 -3 units] of the momentum observables fi3H\ ^3{\) at y/s = 500 GeV for 
three different channels: t + i —> a(q+) + c(q_) + X ; given the Set #2 of SUSY parameters 
(UW). The left column includes only the t (W)EDFF corrections and the right one comes 
from the complete one-loop cross section for e + e~ — > ti. 



CPT 
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b- 
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i- 


- b 
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-t 


even 


Ax 


-0.036 


0.242 
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-0.202 


0.068 


-0.467 


odd 
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-0.180 


-0.204 
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even 


T 33 


-0.006 
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0.021 


-0.140 


-0.037 


0.248 


odd 


Q33 


0.486 


0.542 


-0.335 


-0.374 


-1.274 


-1.420 



over the phase space of the final state particles, 

(0) ac = ~ [<0> aa + (OU\ = J- / [da a - c + dtTes] O , (45) 

2 2cx ac J 

where both the process (ac) and its CP conjugate (ca) are included and 

a ac = Jda a5 = Jda c - a , (46) 

in full analogy with Eqs. (|37],[3^). The differential cross section for t-pair production and 
decay is evaluated for every channel using the narrow width approximation. 

In Table f| the ratio r is shown for three different decay channels at y/s = 500 GeV and 
some realistic CP-odd observables involving the momenta of the decay products analyzing t 
and t polarizations in the laboratory frame. The leptonic decay channels are the best t spin 
analyzers but the number of leptonic events is also smaller. The Reference Set #2 has been 
chosen. As expected, the dipole contributions (left columns) to the CPT even observables 
are very small for this choice of SUSY parameters but the actual expectation values (right 
columns) are larger. 

The previous results were obtained for unpolarized electron and positron beams. Let P± 
be the degree of longitudinal polarization of the initial e . The differential cross section 
reads now 

da = \ [(1 + P+)(l + P_) da R + (1 - P+)(l - P-) da L ] , (47) 

where o~r/l corresponds to the cross section for electrons and positrons with equal right /left- 
handed helicity. Chirality conservation suppresses opposite helicities. Table |5] summarizes 
some extreme cases. If both beams are fully polarized, P + = P = ±1, the ratio r is the same 
as for (P± = 0, P T = ±1), respectively (Table ||), but the cross sections are twice as much 
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Table 5: The same as in Table^ assuming the complete one-loop cross section for e + e — > tt 
and longitudinal polarizations for one of the e beams, P± . 



p± = 


--0, P T 


= -l; o- t f 


= 0.707 pb 


P± = 


-- 0, P T = 


+1; ctt = 


= 0.355 pb 





b-b 


t-b 


l-l 
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b-b 


l-b 


l-l 


At 


0.359 


-0.298 


-0.667 


At 


0.016 


-0.009 


-0.053 


A 2 


0.397 


-0.254 


-0.952 


A 2 


0.134 


-0.099 


-0.468 


T33 


0.065 


-0.214 


0.383 


T33 


-0.003 


0.006 


-0.020 


Q33 


0.708 


-0.488 


-1.848 


Q33 


0.210 


-0.145 


-0.555 



(f!7|), which results in a higher statistical significance of the CP signal. From comparison of 
Table ^ (right columns) with Table [5] is clear that left-handed polarized beams enhance the 
sensitivity to CP-violating effects. 

Concerning the experimental reach of this analysis, the statistical significance of the 
signal of CP violation is given by S = \r\y/N with iV = e£a t tBR(t — > a)BR(£ — > c) where 
e is the detection efficiency and £ the integrated luminosity of the collider. The branching 
ratios of the t decays are BR ~ 1 for the b channel and BR ~ 0.22 for the leptonic channels 
(£ = e, /i). At y/s = 500 GeV the total cross section for t-pair production is a t i — 0.5 pb and 
the NLC integrated luminosity £ ~ 50 fb _1 |21|. Assuming a perfect detection efficiency, one 
gets y/~N ~ 160, 75, 35 for the channels 6 — 6, I — 6, i — I, respectively. With these statistics, 



values of \r\ ~ 10 2 would be necessary to achieve a 1 s.d. effect, which does not seem to be 
at hand in the context of the MSSM, even for polarized beams, as Tables |] and |5] show. 



7 Summary and conclusions 

The one-loop expressions of the dipole form factors of fermions in terms of arbitrary complex 
couplings in a general renormalizable theory for the 't Hooft-Feynman gauge have been given. 
The CP-violating (-conserving) dipole form factors depend explicitly on the imaginary (real) 
part of combinations of the couplings. 

As an application, the electric and the weak-electric dipole form factors of the t quark 
(for s > 4m 2 ) have been evaluated for the MSSM with preserved R-parity and non-universal 
soft-breaking terms. A version with complex couplings of the MSSM has been implemented 
in order to get one-loop CP-violating effects. The results depend on three physical CP 
phases. The supersymmetric parameter space has been scanned in search for the largest 
contributions to the dipoles. The t dipole form factors are larger in the low tan f3 scenario. 
The values depend strongly on the interplay between the energy at which they are evaluated 
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and the set of supersymmetric parameters used as inputs. The real and imaginary parts can 
reach values of similar size: at y^i = 500 GeV, the t EDFF and WEDFF are of O(10~ 19 ) ecm. 

Away from the Z peak both the electromagnetic and the weak dipole form factors are 
equally relevant but not enough to parameterize all the physical effects (in particular CP- 
violation). The case of i-pair production in high energy e + e~ colliders has been considered to 
illustrate this fact. Taking several CP-odd spin-observables based on the t and i polarization 
vectors, the different contributions from vertex and box corrections have been evaluated. 
There is no one loop contribution from the SM sector. It has been shown that, for the set of 
supersymmetric parameters that provides sizeable values of the t CP-violating dipoles, the 
SUSY box contributions happen to contribute with opposite sign and in a similar magnitude, 
yielding altogether a much smaller CP-violating observable effect. Another configuration has 
been shown for which the dipoles are smaller but the combined effect is larger. The same 
analysis has been performed using instead realistic observables based on the momenta of 
t and i decay products, with similar results. The channels in which bb, bl and il act as 
spin analyzers have been used under the assumption of standard CP-conserving decays and 
using the narrow width approximation. As expected, the leptons are the best t spin analyzers 
yielding the maximal values for the signal of CP-violation, r ~ 0.5 x 10~ 3 . Nevertheless the 
statistics for such events is smaller. 
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A Decomposition of the one-loop 3— point integrals 

The systematic way of dealing with one-loop integrals consists of reducing the tensor integrals 
to scalar ones. We employ the following set of 3-point tensor integrals (see Fig. ^|): 



1071" f 

C {0 , ^ } { Pl ,p 2 ,M 1 ,M 2 ,M 3 ) = — /i 4 " D / 

% J 



d D i {l, i„ ijX} 

(2tt) d V 1 V 2 V 3 



(48) 



with 



{I - Pl f - Ml + xe 
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Figure 8: Momentum convention for 3-point integrals. 



£> 2 = (I- P2) 2 - Ml + ie , 

T?3 = I 2 - Ml + ie . (49) 
and the orthogonal reduction (k± = p\ ±p 2 ) |57|1 - 

Cp = C-y k+p + C 1 k_ fl , (50) 

} + C 2 fiV • ( 51 ) 

This decomposition is very convenient because the integral ( f4~8| ) is invariant under the 
combined replacements p\ <-> p 2 and Mi <-> M 2 , and we are dealing with equal mass final- 
state particles p\ = p\ = m 2 . Therefore Cf and C 2 ~ are antisymmetric under these replace- 
ments while the other scalar integrals are symmetric, and the case M\ = M 2 automatically 
yields Cf = C 2 + ~ = 0. 



For comparison we also list the tensor decomposition defined in [58 



% m, M (Pi,P2,M 3 ,M 1 ,M 2 ) = — // y — — ___ , (52) 

with 

V x = I 2 - Ml + ie , 

V 2 = {I + pi) 2 - Ml + ie , 

£>3 = (/+p 2 ) 2 -M 2 2 + k . (53) 

and the tensor decomposition 

^ - Pi M CfWcf>, (54) 



(55) 



These scalar integrals are related to the ones obtained by orthogonal reduction in the fol- 
lowing way: 

C[ D) = -Ct-Ci, (56) 
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d D) = -c+ + cr, (57) 

Cff = Cl (58) 

C[? = Ct + C 2 +2Ct, (59) 

Cg ] = C+ + C 2 -2Ct, (60) 

C[? = C+-CZ, (61) 



with the arguments 

C (D) = C ^)( pi)P2) M3j Ml) Ma)? qj = _ p2? Mi; M2) Ms 



B Conventions for fields and couplings in the SM and 
the MSSM 

B.l The conventions 

We use the conventions of Ref. ||59||. The covariant derivative acting on a SU(2)^ weak 



D, = d, + ig-W; + ig'-B, , (62) 



doublet field with hypercharge Y is given by 

r a Y 
—W a u + ig'- 
2 ^ 2 

where r a are the usual Pauli matrices and the electric charge operator is Qf = /| + Y/2, 
with l[ = r 3 /2. The Z and photon fields are defined by 

Z M = cos 6 W - B^ sin 9 W , (63) 
= W^smOw + B^cosew . (64) 



The charged weak boson field is = (W^ =F iW^) 

In the Standard Model (SM), there is only one Higgs doublet H with hypercharge Y = 1. 
After spontaneous symmetry breaking (SSB), the physical Higgs field H° and the would- 
be-Goldstone bosons x an d ± are given by 

H= U + ^ + w)- (65) 

In the Minimal Supersymmetric Standard Model (MSSM) there are the same matter and 
gauge field multiplets as in the SM, each supplemented by superpartner fields to make up 
complete supersymmetry multiplets. The matter and gauge fields have the same quantum 
number assignments as in the SM. But in the MSSM there are two Higgs doublets with oppo- 
site hypercharges. Each forms a chiral supersymmetry multiplet and an SU(2) doublet. The 
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Lagrangian we use is defined in Ref. |B0]| , especially in Eq. (1.45) and Eq. (1.54). We do not 
impose any reality constraint onto the parameters except for the reality of the Lagrangian. 

Spontaneous breakdown of the SU(2)xU(l) gauge symmetry leads to the existence of 
five physical Higgs particles: two CP-even Higgs bosons h and H, a CP-odd or pseudoscalar 
Higgs boson A, and two charged Higgs particles H ± . They are grouped in two doublets, 
Hi = ir 2 <&l and H 2 = <E>2, with opposite hypercharges (Y = =)=1), where 



$1 



$2 



(66) 



After SSB these are expressed in terms of the physical fields h, H, A, and the would- 
be-Goldstone bosons G°, G^ by 



-H + smf3 + G + cos(3 , 

—j= \v\ + [(—h sin a + H cos a) + \{—A sin (3 + G° cos 
H + cos (3 + G + sin f3 , 

{^2 + [(h cos a + H sin a) + i(A cos (3 + G° sin (3)] | 



(67) 
(68) 
(69) 
(70) 



Besides the four masses, two additional parameters are needed to describe the Higgs sector 
at tree- level: tan/3 = 1)2/1)1, the ratio of the two vacuum expectation values, and a mixing 
angle a in the CP-even sector. However, only two of these parameters are independent. 
Using Ma and tan (3 as input parameters, the masses and the mixing angle a in the H, h 
sector read 



-{M\ + M 2 Z + 6) 



it 



\ 



M\M 2 Z cos 2 2(3 + e{M\ sin 2 f3 + M| cos 2 (3) 



(Mi + Mf + e) 2 



M H ± 
tan 2a 



1 + 



3l 

M\ 



1/2 



tan 2(3- 



M| + M| 



M| — Mf + e/ cos 2/3 
which include the leading radiative correction 



7T 

-- < a < 



3G F / m| N 

v^tt 2 sin 2 /? ° g I + m 2 



(71) 

(72) 
(73) 



(74) 



Here Gp is the Fermi constant and rrig is the common mass scale for the squarks. 

The mass terms for the neutral gauginos and Higgsinos originate from the bilinear Higgs 
part of the superpotential (the \i term) and the soft-SUSY-breaking gaugino mass terms 
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with masses M 2 and M\ for the SU(2) and U(l) gauginos A a and A', respectively. The 
gaugino mass parameters are constrained by the GUT relations 



5 ex 
Mi = - tan 2 6 W M 2 , M 3 = — s 2 w M 2 . 
a 



(75) 



Mixing terms arise from the minimal coupling terms between Higgs, Higgsino and gaugino 
fields, where the Higgs fields have been replaced by their vacuum expectation values. In 
terms of two-component spinors the mass terms add up to 



C£ = -i(-iA',-iA 3 ,<,^ 2 )F(-iA' ) -iA 3 ,<,^ 2 ) T + / i .c. 



(76) 



with the symmetric mass matrix 

/ 

Y = 



Mi 


M z s w cos (3 
\ M z sw sin j3 



. \ 

M 2 

M z c w cos (3 
—M z cw sin (3 —ji J 



(77) 



We define the unitary diagonalization matrix iV and the matrix N', which is often more 
convenient, by the equations 



N*YN~ X = N dmg , N' = N 



( 


c w 


-sw 





\ 




s w 


c w 
















1 





{ 











1 J 



(78) 



In general the \x parameter may be complex and then the neutralino masses depend also on 



the phase of /i. The neutralino mass eigenstates are four Majorana spinors given by 



X] = (P L N jk + P R N* k )*° k , 
-i/ 

iA' 



•3* 
-iA 3N 

iA 3 . 



A, 



(79) 
(80) 



The mass terms for the charged gaugino and Higgsino have the same origin and, in terms 
of two-component spinors, read 



rx 



X = 



M 2 

M w y/2 cos 13 




(81) 



with the abbreviations 



^ + = (-iA + ,^+ 2 ), V7 = (-iA-^H 1 ) 



(82) 
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We now define unitary diagonalization matrices U, V by the equation 

WXV- 1 = M diag . (83) 
The masses of the charginos are given by 

m 2 ± =l{ M 2 2 + \fi\ 2 + 2M^ T \{M 2 2 - \fi\ 2 ) 2 + 4M^ cos 2 2(3 (84) 
+4M^(M 2 2 + |/i| 2 + 2M 2 Re(/i) sin2/3)] 1/2 }. (85) 

The chargino mass eigenstates are two Dirac spinors \j given by 

Xt = (PLV jk + P R U* k )y k: (86) 

We abbreviate the charge conjugated fields as 

Xi = Xt° ■ (88) 

The mass terms for the scalar quarks originate from the Yukawa couplings to the Higgs 
fields (which yield the corresponding quark masses) and the soft-SUSY-breaking squark 
mass terms and squark-Higgs interactions parameterized by rrig L/R and A q , respectively. 
Moreover there are mass and mixing terms from auxiliary field terms involving one or two 
Higgs bosons and two squark fields. The resulting mass matrix for the two squarks of flavor 
q is: 



, , i m q T 2 + m 2 m 9 T *p>m a 

M\= q q q \\ \ , (89) 

m LR m q m R + m " 



where 



m 
m 



f = ml +C os2(3(li-Q f s 2 w )M 2 z , (90) 
I 2 = m 2 + cos 2(3 {Q f s 2 w )M 2 z (91) 



and 



m\ R = A q - //*{cot (3, tan (3} , (92) 



for {up, down}- type squarks, respectively. 

This hermitian mass matrix is diagonalized by a unitary matrix S q , so we can write the 
squark mass eigenstates of flavor q as 

qi = Sn qi + Sf 2 q R , 

h = S q 21 q L + S q 22 q R . (93) 
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In general, va\ R = \mL R \e liP9 is a complex number and S g may be written as 

sq= ( eW»cos0 8 e-^ainflj \ 
" ^ - e W2 s i n 0. e -W2 cos ^ y 

The mass terms of the sleptons arise in the same way as the squark masses. The main 
difference appears for the sneutrinos: there is only one sneutrino for every generation, fy, 
and hence there is no mixing. Moreover, we cannot add trilinear soft-breaking terms to shift 
the masses of the sneutrinos, whose value is given just by m~ = m| + 1/2 cos 2/3 M|. (See 
Eq. ©). 

In general, the physical sfermion masses are given by 



I 1 2 2 / 2 2 1 

m /i, 2 = m ) + o 1 ( m £ + m i ) =F V ( m R ~ m L ) 2 + ^)\m{ R \ 2 \ , (95) 



independent of the phase of 



B.2 Physical Phases in the MSSM with complex couplings 

As we do not constrain the parameters of the MSSM to be real, the following parameters may 
take complex values: the Yukawa couplings, fi, the soft parameters mf 2 , M\, M2, M3 and 
the A parameters. But not all combinations of phases in these parameters lead to different 
physical results because several phases can be absorbed by redefinitions of the fields. We 
now describe a procedure to eliminate the unphysical phases. We assume the GUT relation 
between the Mj, so they have one common phase. The only remaining phases will be chosen 
to be those of /i, the A parameters and, for three generations, one phase for all the Yukawa 
couplings (the <5ckm)- 

Analogous to the Standard Model case, the Yukawa couplings can be changed by redef- 
initions of the quark superfields in such a way that there remains only one phase for three 
generations and only real couplings for less than three generations. 

Table 6: The charges rii for three U(l) symmetries that leave £mssm invariant. 



U(l) 




A 


m\ 2 




Hx 


H 2 


QU 


QD 


LE 


e 


PQ 
Ri 



-1 



-1 


-1 



-1 
1 


1/2 



1/2 



-1/2 
1 


-1/2 
1 


-1/2 
1 




-1/2 



Following |12| we consider two U(l) transformations PQ and R\ that do not only trans- 
form the fields but also the parameters of the MSSM. In table ^| the charges ni of the various 



32 



quantities are displayed with which the Lagrangian £mssm is invariant under the combined 
multiplications with e tan \ The first transformation is a Peccei-Quinn symmetry and R\ 
is an R symmetry that also transforms the 9 variables appearing in the arguments of the 
superfields in £mssm- 

Since £mssm is invariant under these transformations, so are the physical predictions of 
the MSSM. If the parameters Mi, A, m\ 2 and \i have the phases <pm, <Pa, fm^ an d </V> we 
first apply R± with the angle pm, then PQ with the angle (p m ^ 2 and obtain for an arbitrary 
observable: 



= <r(\p\,\ A l\ M il\ m i2l <P» + 9m, <Pa-<Pm, 0, (96) 

= °"(I/^U^UM|, |m? 2 |, <Pp + <pM-<Pm* 2 , (PA-fM, 0, 0). 

So the physical predictions only depend on the absolute values of the parameters and the 
phases 

A = arg(AM*), B ee arg(/Mm?*) . (97) 

One can replace these phases by another set where only the \x and the A parameters are 
complex. In our choice, the phase of A is traded for the phase tp? of the off-diagonal term 
in the corresponding sfermion mixing matrix: m LR = Af — /x*{cot, tan}/3 (|92|). Relaxing 
universality for the soft-breaking terms, every Af contains a different CP-violating phase. 

B.3 Vertex factors 

We employ the notation and conventions of the previous sections. The list of generic vertices 
is shown in Fig. ||. 

B.3.1 Couplings in the SM 

Vertex 1: Coupling of electroweak gauge bosons and fermions: iey^iV — Ay^). 

iff - V =-Q f , A = 0, (98) 
v f A _ a f 

1-SwCw IswCw 

1 

2\/2sw 

with v f ee (// - 2s 2 w Q f ), a f = //. 
Vertex 2: Coupling of gluons and quarks: — ieg s 7^T a 



Zff'- V =-^-, A — —777^— , (99) 
W ± f'f: V =A = --±—, (100) 
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Figure 9: Generic Vertices. 



Vertex 3: Coupling of two fermions and one Higgs boson (Tiff): ie(S — -P75) 



with fif 



H°ff 

xff 

TO/ 



S = -n f , P = 0, 

S = 0, P = -2Uln f , 

S = y/2li\pL f - n r \ , P = -V2li[fi f + n f ,\ , 



2s w M w 2s w c w M z 
Vertex 4: Coupling of one Higgs boson and two gauge bosons: ieKg^ v . 



H°ZZ : 


K = 


M z 






SwCw 


H°W ± W T : 


K = 


M w 








(j) ± W T -f : 


K = 


M w , 


(j) ± W T Z : 


K = 


-M z s w 



The rest of the couplings are K — 0. 

Vertex 5: Coupling of one gauge boson and two Higgs bosons (V<pi(p2)- 
ieG(p 1 +p 2 ) tl . 



(101) 
(102) 
(103) 



(104) 

(105) 

(106) 
(107) 



70 A <T A : G= -A, 



(108) 
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Z(f) X (f)~ X : G 

Z X H°: G 

W x <p' x H G : G 

W x (j)~ x x- G 



-A 



cos 29 w 

2sy/Cw 



A 



2sw 
i 



2s 



w 



(109) 
(110) 
(111) 
(112) 



Interchanging the two Higgs bosons causes the coupling constant to change sign. 
Vertex 6: Coupling of three gauge bosons (outgoing momenta): 

ieJ{ 9llu (k x - k°) p + g up (k- x - k x ), + g Pfl (k° - k~ x ) v } . 



1 w x w- x 
zw x w- x 



J 
J 



-A, 



-A— 



(113) 
(114) 



B.3.2 Couplings in the MSSM 

Vertex 1: Coupling of two neutralinos or charginos and one gauge boson: icy^V — A75). 

The fermion flow direction in our Feynman graphs is fixed by the outgoing fermion 
antifermion pair. The relation between the gauge boson vertex factors for the two dif- 
ferent fermion flow directions is obtained by substituting spinors by charge conjugated 
spinors in the interaction term: 



V^ + Y(Pl9l + Pr9r)X + = Y(~Pl9r ~ Pr9l)x 



(115) 



Taking into account also the symmetry factor S = 2 for the neutralino coupling and 
S — 1 for the charginos, the boson vertex factors V = (g L + g R )/2 and A = (g L — g R )/2 
are given by: 



,7-0 ~o 

ZXjXk ■ 9l 



ZXk x\ 



ZXj X k 



9r 
9l 
9r 
9l 
9r 



2s\yC\V 
1 

2s\yCw 
1 

SwCw 

1 

SwCw 

1 

SwCw 

1 

SwCw 



(N' k4 N'* 4 - N^N") , 
(KlN' j3 - N'C 4 N> 4 ) , 



' w 

,2 



V k \V j2 + c w V k \V n 
2 



2 - s w j U k2 U* 2 + c w U kx U] x 

1 _ s 2 \ TT TT* 1 „2 

2 

4^ + iKl^l 



s z w \u k2 u; 2 + c^u kl u; 1 
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(116) 
(117) 
(118) 
(119) 
(120) 
(121) 



and 



IXk Xj ■ 9l=9r= -S jk , 
IxJXk 9l=9r = S jk . 



(122) 
(123) 



Vertex 2: There is no genuine supersymmetric vertex of this kind. 

Vertex 3: There are three couplings of the kind ie(5' — P75) : 

(Some of the couplings are more easily written in terms of g^R = S ± P). 
• Coupling of one neutralino or chargino to fermions and scalar fermions. 
The couplings of neutralinos and charginos to quarks and scalar quarks are given by 

C>xm = etfZ(P L g L + P R g R )q + eq(P R g* L + P L g R )xq • (124) 

down-type quarks (i = 1,3,5): 



QtkXjQi ■ 9l 



= -V2 



9r 



9l 
9r 



up-type quarks (i = 2,4,6): 



vtkXjQi ■ 9l 



9r 



ii-i),kxUi : 9l 



9r 



QiN'A + (P 3 - Q lS 2 w ) 



1 



SwCw 



Nf 2 ) SI, 



AT'* Qi 

2M w s w cos j3 k2 



-V2 



-(Q iN '+{-Q iS 2 w ) 



Sy/Cw 



N' J2 S\ 2 



rn 



2M w s w cos/3 

1 



-nss 1 + 



m 



1 

sw 



V2M W cos (3 



Qi+l -y* 

V2M w smp j2 k2 



-V2 



QtN* + (P 3 - Q t s 2 w ) N' j2 j u kl 

s w c w 



2M w s w sin/5 



|\f" Q* 



= -V2 



QiN' + (-Q iS 



WJ 



SwCw 



N' J2 ) S\ 2 



2M w s w sin/3 

_L _/7* 4. 
sw 

1 m„ 

sw 



Jv j4°fcl 



V2M W sin 



y/2Mw cos f3 j2 k2 . 



V^ 1 



(125) 



(126) 
(127) 
(128) 



(129) 



(130) 
(131) 
(132) 
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The couplings of neutralinos and charginos to leptons and scalar leptons are given 
analogously, performing the following substitutions: 



1, 3, 5 : qn — > l\ , S* k — > S l : k , Qi — — 1 , Jg — — — , m q 



Qi2 — > h , 

2,4,6 : q a -> i>, , Sj fe -> <J jfc , = , P 3 = - , m 9i = , 
§i2 does not exist . 



• Coupling of one gluino to a quark and a scalar quark. 
The interaction between gluinos, quarks and squarks is described by the terms 



C = $e~g (P L g L + P R g R )-q + eq(P R gl + P L g* R )-~9 a q i 
yielding the vertex factors 

ie(P L g L + P R g R ) — , ie(P R g* L + P L g* R ) — . 
In our calculations the Gell-Mann matrices appear only in the combination 



AB 



C 2 (F)Sab = 3^b- 



The couplings are 



Qi9Q e 9L = -V2g s S^ , 



egn. 



+V2g s S,, 



Coupling of two fermions and one Higgs boson (Tiff)- 



Huu 
Hdd 

huu 

hdd 
Auu 
Add 

H + ud: S 



S = —fJ> u sin a I sin (3 . 
S = —/id cos a I cos (3 
S = — /i« cos a/ sin (3 . 
S = fid sin a I cos (3 , 

S = , 

S = , 



P = 0, 
P = 0, 



p 
p 
p 
p 



o, 
o , 

— ifi u cot (3 
—i/id tan (3 



1 

72 



(fj, u cot (3 + fid tan f3) 



G°ff : 
G ± /7 : 



P = ^(/i u cot/3 - tan/5) 
S = , 
S , = v^/ 3 / ^/-A*/'] , 



P 
P 



-21/l/i/ , 
-V2l(\n f + /i// 



(133) 



(134) 



(135) 



(136) 
(137) 



(138) 
(139) 
(140) 
(141) 
(142) 
(143) 

(144) 
(145) 
(146) 



with fj,f = rrif /2swM w = m,f /2swCwM Z - For the vertices corresponding to the her- 
mitian conjugated fields we have to replace (S, P) by (S*, —P*). 



37 



Vertex 4: Coupling of one Higgs boson and two gauge bosons: leKg^. 
The couplings of a neutral Higgs to two Z bosons are: 

hZZ: K= — sin(/3 - a) , (147) 

s w c w 

My 

HZZ : K = — cos(/3 - a) . (148) 

swcw 

And the only nonzero couplings of a charged Higgs to neutral gauge and W bosons 
are: 

G ± W T -f : K = M w , (149) 
G ± W T Z : K = -M z s w . (150) 

The rest of the couplings are K — 0. 

Vertex 5: There are two types of couplings Z0102 of the kind ieG(pi +P2)n ■ 
• Coupling of two Higgs bosons and one neutral gauge boson. 



(151) 

(152) 

(153) 

(154) 

(155) 

(156) 
(157) 



ZAH : 


G = 


isin(/5 — a) 

2s\yCw 


ZAh : 


G = 


i cos(/3 — a) 


ZG°H : 


G = 


icos(/3 — a) 

2s\yCw 


ZG°h : 


G = 


i sin(/3 — a) 

2syyC\y 


ZH X H- X : 


G = 


^cos2^ 

2s\yCw 


ZG X G~ X : 


G = 


^cos 29 w 

2-S\yC\y 


-fH x H- x : 


G = 


-A, 


~lG X G- X : 


G = 


-A . 



(158) 

Interchanging the 0i and 2 causes the coupling constant to change sign. 
• Coupling of two scalar fermions and one neutral gauge boson. 

lQijQi,k '■ G = -Qi$jk , (159) 
Z%jQi,k '■ G = (-^3 ~~ Qi s w)Sj*Ski ~ Qi s w^j*2^k2 ■ (160) 

Vertex 6: There is no genuine supersymmetric vertex of this kind. 
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